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We demonstrate a possibility to create a new state of ultracold atoms which we call a Bose-Einstein- 
Young condensate. Atoms are supposed to be in different hyperfine state of the same isotope. The 
wave function of such a state, although totally symmetric with respect to simultaneous permutation 
of co-ordinates and spins of any pair of atoms, has more complicated structure than a simple product 
of totally symmetric co-ordinate and spin parts. Its properties with respect to permutations of only 
co-ordinates or only spins are characterized by a particular Young diagram, a symbol denoting an 
irreducible representation of the permutation group. 
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Creation of binary mixtures of Bose-Einstein conden- 
sates (BECs) of ensembles of atoms of the same isotope 
in different hyperfine states [1,2] and spinor condensates 
in optical traps [3,4] is one of the remarkable advances 
in experiments with BECs. The results of these experi- 
ments are in good agreement with the existing theories 
of a multicomponent BEC [5,6]. Even if one deals with 
a spinor BEC, where an external magnetic field gives 
rise to magnetization and, hence, non-trivial spin prop- 
erties, including condensate fragmentation [6], the trans- 
lational motion of atoms at sufficiently low temperature 
(disregarding small quantum depletion) is essentially the 
ground state of a trap that can be determined by energy 
minimization. 

However, a question arises, is such a ground state al- 
ways a final state reached at the end of a process of cool- 
ing of atoms. In the present Rapid Communication, we 
demonstrate that the answer is negative. It is possible 
to find various states those should appear as final states 
of a sympathetic cooling process and are not coupled by 
interatomic interactions to the ground state determined 
in Refs. [5,6]. 

Due to the Bose-Einstein statistics, the wave function 
of the atoms must be invariant with respect to permuta- 
tion of both the co-ordinates (rj) and spins (<x,-) of any 
pair of atoms. In general case, it can be constructed in 
different ways, namely, 
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Here {A} denotes an irreducible representation of the 
symmetric group of the order n (i.e., the group of permu- 
tations of n objects), n is the number of identical atoms 
in the system, and is the dimension of the represen- 
tation. The theory of symmetric group representations 
and its applications in physics are described, e.g., in clas- 
sical textbooks by Wigner [7] and by Lyubarskii [8]. It 
is worth to note that theory of symmetric group repre- 
sentation is widely used in theory of atomic spectra [9] 



as well as of rotational and vibrational spectra of sym- 
metric polyatomic molecules containing identical nuclei. 
Here we remind briefly the necessary subjects from the 
theory of symmetric group representations. 

Any irreducible representation of the symmetric group 
of order n is characterized by a partition {A} = 
{Ai, A2, ... , A r }, i.e., representation of the number n as 
a sum of integer positive summands arranged in the non- 
increasing order: 

n = Ai + A 2 + ... + A r , 



Ai > A 2 > ... > A r > 0. 

Graphically the partition {A} can be represented by a 
so-called Young diagram consisting of n square boxes ar- 
ranged in rows, the number of boxes in z-th row is equal to 
Ai. To get a basis for this representation, one has to asso- 
ciate each co-ordinate or spin variable to one of the boxes. 
Then it is necessary to apply to a co-ordinate or, re- 
spectively, spin wave function a certain symmetrization- 
antisymmctrization procedure that involves particular 
permutations of the function arguments, according to 
their distribution between the boxes, and summation of 
the obtained terms with the sign of a term depending on a 
given permutation operator. This procedure is explicitly 
described in various texbooks, e.g., in Ref. [8]. The num- 
ber of linearly independent functions obtained in such a 
way starting from different distribution of variables be- 
tween the boxes of the Young diagram is the dimension 
of the representation. It can be calculated using the fol- 
lowing formula [9]: 
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where ij = Ai +r — i, i = 1, 2, ... , r. The sum of squares 
of dimensions of all the irreducible representations equals 
to the number of operations in the group: 
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A basis of a given representation can be orthogonal- 
ized and normalized (by means of linear transformation) 
with respect to the usual quantum mechanical definition 
of the scalar product, implying integration over coordi- 
nates and summation over spin variables. Thus we get 
the two sets of functions <j>^ and Xk^*> depending on 
atomic co-ordinates and spins, respectively. If we denote 
the operators of permutations of the co-ordinates and 
spins of the z-th and /-th atoms by P(i, i') and Q(i, i'), 
respectively, the transformation law can be written as 
follows: 

P(M0<^ A} = £T fem (M')<4 A} , (4) 
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Note that the same matrix Tk m {i, i') appears in 
Eqs.(4, 5). Since the scalar product of quantum mechan- 
ical wave functions is invariant under the permutation 
operations, and the chosen basis is orthonormal, this 
matrix is unitary, J2k T krn(i, i')T^ m ,(i, i') = 8 mm , . Ob- 
viously, the functions xj^ complex conjugate to x^ * 
also form a basis of an irreducible representation (a con- 
jugate representation), and the transformation matri- 
ces are Tj* m (i, i'). Hence, the wave function ^ given 
by Eq.(l) remains unchanged after application of opera- 
tor P(i, i')Q(i, i') permuting both the co-ordinates and 
spins of the pair of atoms. Also it is proved in the group 
theory that there is no way to construct the wave func- 
tion having the Bose-Einstein symmetry property be- 
sides that given by Eq.(l). However, the explicit form 
for the bosonic wave function in a case of an arbitrary 
{A} is extremely lengthy. As an example, we write the 
formula for three particles, {A} = {2, 1}, a, b being the 
translational motion states and 1, 2 being the spin states: 

* = 371 ^ a ( ri ^ fc ( r2 ) - V>a(r 2 )^&(ri)]Va(r3) 
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Consider an atomic ensemble composed of N atoms 
in the state |1) and N atoms in the state |2). The 
states |1) and |2) can be identified with the states 
\F = 1, nip = —1) and \F = 2, nip = 1), respectively, of 
the ground state of the same bosonic isotope, e.g., 23 Na 
or 87 Rb. These states have equal magnetic momenta and, 
hence, experience the same trapping potential. For the 
sake of simplicity, we assume N to be small enough to 



allow us to use a single-particle classification of trans- 
lational motion states in the trap. Energy levels are 
essentially the levels of a single atom in the trap, the 
interatomic interaction provides only small correction to 
them. This assumption is not necessary for validity of the 
subsequent treatment, however, it allows us to present 
our arguments in the most compact and brief manner. 
The states of translation motion will be numbered by j, 
j = 0, 1, 2, in the non-decreasing order in energy, so 
that means the ground state of the trap. The tilde 
symbol is used to distinguish these states from the spin 
states |1) and |2). We use also notation n = 2N for the 
total number of atoms in the trap. 

In the standard theory of a multicomponent BEC, even 
when dealing with the elementary excitations of BEC, the 
wave function is always assumed to correspond to the 
simplest possible Young diagram characterizing both the 
co-ordinate and spin parts. This diagram has only one 
row consisting of n boxes, where n is the total number of 
atoms in the system; the dimension of the corresponding 
(totally symmetric) irreducible representation of the per- 
mutation group is 1, so that usually one writes the wave 
function as a product of co-ordinate and spin parts, 

*(ri, ... , r„; 01, ... , a n ) = </> W (ri, ... , r„) 

*X {n} (<Ti,-,<r n ), (7) 

and the functions in the right-hand-side of this equation 
are invariant with respect to permutations of only co- 
ordinates or only spins. Now we show a way leading 
to creation of complicated states corresponding to the 
general form of Eq.(l). 

Initially, before the cooling process begins, atoms are 
thermally distributed. Mean occupation number for each 
translational state is much less than unity, so that we can 
assume that all the trap states occupied by the atoms are 
different. 

The system of thermal atoms is not in a pure state, 
but rather is described by a statistical mixture of states 
characterized by all the possible Young diagrams. The 
statistical weight of each Young diagram in this mixture 
is directly proportional to the number of different basis 
sets of functions those can be associated with the dia- 
gram. The basis sets arc different in the sense that in 
any subspace of the total Hilbert space spanned by each 
basis there are Hilbert vectors not belonging to any of the 
remaining subspaces spanned by other basis sets. This 
is rather a subtle point. When we discussed earlier a 
question of dimensionality of an irreducible representa- 
tion of the permutation group, we dealt with permuta- 
tion of only function arguments within a Young diagram 
and assumed implicitly that each single-particle state j 
was associated with a fixed box of the diagram. However, 
to determine a statistical weight of any Young diagram, 
we have to allow also permutation of state indices within 
the diagram. As it was assumed before, the probability 
that a given state of translational motion is occupied by 
more than 1 atom is negligible before the cooling process 
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starts. And the number of non-trivial ways to assign n 
different state indices to a Young diagram consisting of n 
boxes equals exactly to the number of non-trivial distri- 
bution of n co-ordinate variables between the boxes of the 
same diagram, i.e., to s^ x \ So that for a given {A} there 
are ways to construct different s^-dimensional ba- 
sis sets of co-ordinate functions. All this sets together 
comprise, according to Eq.(3), n! functions [cf. Eq.(3)]. 

However, not all of the Young diagrams appear in our 
problem. The restriction on possible Young diagrams 
arises from that fact that only two spin states are present 
in this Bose gas mixture. Only diagrams consisting of two 
rows, at maximum, can appear, and {A} = {2N — K, K}, 
K = 0, 1,2, ... , N. If one tries to construct spin wave 
function of N atoms in the state |1) and N atoms in the 
state 1 2) characterized by a Young diagrams with 3 or 
more rows, the result will be identically zero [7,9]. Hence, 
only Young diagrams with not more than two rows will 
apply also to functions of co-ordinates in our problem. 

There is only one non-trivial way to assign a spin state 
label (1 or 2) to each box of the diagram, so that spin 
effects do not increase the statistical weights of diagrams. 

Thus, the statistical weight of the Young diagram 
{2N — K, k} in the initial thermal mixed state of the 
atomic ensemble is, according to Eq.(2) (see also Ref. 
[7]), equals to s i 2N - k ^} = C? N - C^ N \ where = 
n!/[(n — K)\K\] is the binomial coefficient, and K = 
0, 1, ... , N . The maximum of this probability distribu- 
tion locates very close to K = N . Calculations show 
that the average length of the second (shorter) row of 
the Young diagram will be approximately (for N 3> 1) 

(K) & N — y/irN/2. (8) 

In other words, most probably, an observation will result 
in finding a state characterized by a Young diagram with 
two almost equal rows, which is very far from the Young 
diagram with only one row, and the probability to get 
the latter diagram tends to zero when TV increases. 

Then consider a cooling process. As it was said before, 
we assume a sympathetic cooling of our two-component 
atomic ensemble due to interaction with the reservoir of 
atoms of a different isotope kept at a sufficiently low tem- 
perature. 

An interaction operator acting only to atomic co- 
ordinates or only to atomic spins cannot couple a total 
wave function of n identical bosons characterized by the 
Young diagram {A} to a wave function characterized by 
{A'} 7^ {A}. Even an interatomic potential depending on 
an atomic spin orientation cannot cause such a transi- 
tion during an elastic collision. Only if both the spin and 
translational motion states of an atom change, a tran- 
sition between states characterized by different Young 
diagrams becomes possible. But such a process is, in 
fact, a spin-flip collision leading to inelastic losses but 
not facilitating the cooling process. The cooling can be 
accomplished successfully if the inelastic processes are 
inefficient on the cooling time scale. 



Thus we can see that during the sympathetic cooling 
process the type of the wave function symmetry with re- 
spect to permutations of co-ordinates only is conserved. 
This means that at the final stage of the cooling process 
an observer also will find mostly states characterized by 
Young diagrams with two rows of almost equal length. 
But the translational motion wave functions correspond- 
ing to such diagrams cannot be constructed using only 
ground state of the trap as a one-particle state appear- 
ing n times. If one tries to construct a co-ordinate wave 
function by assigning to every box of a two-row diagram 
the same state, one will end up with identical zero. 

Therefore, the lowest energy state corresponding to the 
Young diagram {2N — K, K} corresponds to filling of 
2N — K boxes in the first row of the diagram with the 
ground one-particle state of the trap (|0)) and K boxes in 
the second row of the diagram with the first excited state 
(|l)). Thus the total energy is 2NEq+K(Ei~Eq), where 

Ej is the energy of the j-th one-particle state. Usual BEC 
(corresponding to the Young diagram {2N}) would have 
the energy 2NE^. 

Thus we see that the evaporative cooling process of 
a random mixture of the same isotope in the hyperfinc 
states |1) and |2) with almost 100 % probability results 
in creation of a new state of degenerate atomic Bose gas, 
where almost half of atoms are in the first excited state 
of the trap, and no further cooling can put them into 
the ground state. We propose to call such a macro- 
scopic quantum object a Bose-Einstein- Young conden- 
sate (BEYC). 

We would like to stress once again that BEYC can be 
converted to an ordinary BEC only by processes where 
both the spin and translation motion of atoms is changed. 
But spin-flip collisions being the most obvious example 
of such a process lead to losses and heating of the atomic 
sample, due to hyperfine energy release. Spin-flip colli- 
sions, hence, simply destroy BEYC as well as BEC rather 
than facilitate transitions between them. 

The situation with BEYC and ordinary BEC is similar 
to that of orth- and parahydrogen molecules. Transition 
between these molecular states also requires change in 
the type of a Young diagram ({2} or {1, 1}) and, hence, 
requires a rather long time even for room temperature 
and atmospheric presssure. 

The number of atoms in the first excited trap in the 
BEYC state far exceeds both the quantum and thermal 
(at experimentally reachable temperatures) condensate 
depletion. The first excited state is, in fact, the dipolc 
oscillation state. If the number of atoms in each of the 
two hyperfine states is well defined, the phase of these 
oscillations is indeterminate. However, it must become 
well-defined if the BEYC undergoes a continuous mea- 
surement of atomic positions, in close analogy with for- 
mation of an interference pattern produced by two coun- 
terpropagating condensates each of them was initially in 
the Fock state [12]. 

Now we have to say few words on experimental per- 
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spectives for a BEYC. First of all, we must answer, 
why BEYC has not been obtained up to now. In Ref. 
[1], the evaporative cooling technique was used. How- 
ever, 87 Rb was present in that experiment in the states 
\F = l,m F = -1) and \F = 2, m F = 2). Their mag- 
netic momenta differ by a factor 2. Hence, they experi- 
ence different trapping potentials, and, due to the grav- 
itation, the centers of the potential wells are displaced 
one from another. The displacement is of such a mag- 
nitude that, on the cooling stage, the thermal clouds of 
atoms in different spin states overlap significantly, but 
the two BEC do not. In fact, condensation of rubidium 
atoms takes place in different states for different mp's. 
In such a situation, it is hard to distinguish between the 
states characterized by different Young schemes, since all 
of them have the same energy, N_\Eq _ 1 +N 2 Eq 2 , where 
N mF is the number of atoms with the magnetic quantum 
number rap and Eq mp is the energy of the ground state 
of these atoms in the corresponding trapping potential 
influenced by the gravitation. In Ref. [2] an ordinary 
(single-component) BEC was prepared first. Than a frac- 
tion of the atomic ensemble was transferred to another 
hyperfine state by a two-photon microwave- RF pulse. Al- 
though the magnetic momenta of the two states were the 
same, creation of a BEYC was prevented by that cir- 
cumstance that the transferred kinetic momentum (cor- 
responding to the hyperfine splitting plus Zeeman shift) 
was negligible compared to the inverse trap size. Thus 
the translational state of atoms was not changed, and, 
hence, the wave function retained its symmetry prop- 
erties with respect to permutation of only co-ordinates 
or only spins (i.e., remained a product of totally sym- 
metric co-ordinate and spin parts). In the experiment 
of Ref. [3], the atoms were prepared initially in the or- 
dinary (single-component) BEC state in a optical trap, 
and then the atomic spins were allowed to evolve due 
to interatomic interactions. The BEYC state was not 
reached there since it had a significantly higher energy 
than that stored in the BEC. The most obscure case is 
the case of all-optical formation of a BEC [4] . In this ex- 
periment, an evaporative cooling technique was used, and 
the number of atoms in BEC was more than by order of 
magnitude smaller that the number of atoms just before 
the beginning of the cooling process. This case is be- 
yond the scope of the present Rapid Communication and 
requires straightforward numerical (Monte Carlo) mod- 
eling. Only performing such a modeling, it is possible to 
understand how Young diagrams, in average, change. A 
scenario, when more than 90 % of atoms are lost dur- 
ing the evaporative cooling process, complicated Young 
diagrams do not "survive" , and the Young diagram con- 
sisting of only one row of boxes prevails on the final stage, 
looks quite reasonable. 

Thus, to observe BEYC one needs to repeat the exper- 
iment like that of Ref. [1] (where the atoms in the two dif- 
ferent hyperfine states were prepared independently), but 
with the states \F — I, mp = —1) and \F = 2, to_f = 2) 
those have equal magnetic momenta. The combined mag- 



netic trapping and gravitational potential will be the 
same, no spatial separation of condensed atoms in dif- 
ferent hyperfine state will occur, and difference of BEYC 
from ordinary BEC will become apparent. An alterna- 
tive approach would imply transfer of atoms into another 
hyperfine state not by a microwave- RF pulse as in Ref. [2] 
but by a Raman transition in two crossed laser beams. A 
non-zero angle of crossing will provide a kinetic momen- 
tum transfer needed for a translational motion change. 
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